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Abstract 

In this talk I consider QCD sum rules for the ground state heavy baryons to 
leading order in 1 /mg and at next-to- leading order in as within the context of Heavy 
Quark Symmetry. The analysis is done at a fixed scale /i = 1 GeV. The evolution 
behaviour of the residues is determined by the two-loop anomalous dimension of the 
heavy baryon currents calculated earlier. I compare the diagonal, non-diagonal and 
constituent type QCD sum rules and show that the diagonal sum rules are the most 
reliable one. As central values for the bound state energies I find to(Aq) — tuq ~ 
760 MeV and m(T,Q) — tjiq ~ 940 MeV. The central values for the residues are given 
by F(A Q ) ~ 0.030 GeV 3 and F(Eq) ~ 0.038 GeV 3 . 

1 Talk given at the IVth International Workshop on Progress in Heavy Quark Physics, 
Rostock, September 20-22, 1997 
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1 Introduction 



The correlator of two baryonic currents is the fundamental building block for stating QCD 
sum rules for heavy baryons, i.e. baryons which contain one heavy quark. In this talk I 
present the calculation of first order radiative corrections to QCD sum rules within the 
Heavy Quark Effective Theory (HQET) in the limit of infinite mass for this quark . The 
calculation of first order radiative corrections 0, |3| make sense because of the knowledge 
of the two-loop anomalous dimension || . The interpolating currents for the heavy baryons 
can be written as 

Jb = [(gi) jT CTr(g 2 )^]r'g% fc , (1) 

where the index T means transposition, C is the charge conjugation matrix, r is a matrix 
in flavour space, and i, j and k are color indices. V and V are the light-side and heavy- 
side Dirac structures of the current vertices. For each ground-state baryon there are two 
independent current representations || carrying the same quantum numbers, for the Aq- 
type baryons they are e.g. given by 

J Al = [(qi) iT Cr l5 (q 2 y}Q k e ijk and J Aa = [{q^Cr^q^Q^^. (2) 

2 The two-point correlator 

The two-point correlator can be constructed as 

U 13 (uj =p. v )=ij (0|T{ Ji(x) J i (0)}|0)e ip Vx = T'^-f itr(rf )2tr(rr+)^(w). (3) 

i and j stand for the two independent current representations for the ground state baryons. 
This gives rise to two independent types of correlators, namely the diagonal correlators 
(Ji J\) and (J2J2) and the non-diagonal correlators (J1J2) and (J2J1) as well as linear com- 
binations of these cases. Following the standard QCD sum rule method || , the correlator 
is calculated in the Euclidean region — uo ^1 — 2 GeV including perturbative and non- 
perturbative contributions. In the Euclidean region the non-perturbative contributions are 
expected to form a convergent series. The non-perturbative effects are taken into account 
by employing an operator product expansion (OPE) for the time-ordered product of the 
currents. One then has 

(T{J(x),J(0)}) = Y,C d {x 2 )O d = (4) 

d 

= C (x 2 )O + C 3 (x 2 )0 3 + C 4 (x 2 )0, + C 5 (x 2 )0 5 + ... 

where the Od are vacuum expectation values of local operators whose mass dimensions 
are labelled by their subscripts d. Oq = 1 corrersponds to the so-called perturbative 
term, 3 = (qq) is a quark condensate term, 4 = a s (G 2 ) is a gluon condensate term, 
5 = Qsiq&nvG^q) is a mixed quark-gluon condensate, and so on. The expansion coefficient 
Cd(x 2 ) are the corresponding Wilson coefficients of the operator product expansion. The 
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different diagrams are shown in Fig. 1. 

A straightforward dimensional analysis shows that the operator product expansion of the 
diagonal correlators contain only even-dimensional terms, while the expansion of the non- 
diagonal correlators contain only odd-dimensional terms. This classification is preserved 
when radiative corrections are included, assuming the light quarks to be massless. 



3 Construction of QCD sum rules 

The QCD sum rules can be constructed by taking care of two possible expressions for the 
two-point correlator. On the one hand, the scalar correlator function P{uj) satisfies the 
dispersion relation 

P(u) = / 7 v ' + subtraction, (5) 

Jo w'-w-jO 

where p{uj) = lm(P{uj))/n is the spectral density. On the the phenomenological side of 
view, the two-point correlator is represented by the spectral representation 

U(u) = - V l( ° |Jx|X)|2 + subtraction. (6) 



The scalar correlator funktion can thus be expressed by the residues Fb given by 

1 

71 



J A \A Q ) = F A u, (0|J s |E Q ) =F^u and (0|J£*|E*) = ^=F^u v (7) 



in the form 12 12 

P( u ) = - + y 2^1 + subtraction ( 8 ) 

A — uj — iO ujx — uj — iO ' v ; 

where A = m B — Mq is the ground state energy of the baryon. The main assumption of the 
QCD sum rule method is that the remaining sum can be approximated by the integral of the 
spectral density given by the dispersion relation starting from some threshold energy Eq- 
The combination of the phenomenological and the theoretical identity for the correlator 
function leads to 

7 21 = / y ' + subtraction. 9 

A-u-jO Jo w'-w-jO 

This formula is not useful anyhow, since the spectral density is reliable only for negative 
values of uj, while the integral is to be calculated mainly at uj = A. This region of integration 
can be reached by an extrapolation using higher and higher derivatives when uj goes to 
—00. This extrapolation is expressed by the Borel transformation 

BU T ) = bP(](uj)):= lim ^ '- ^—f(ou), T = fixed. (10) 

J v ' uj,n^oo n \ duj n v ' n v ' 
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This Borel transformation leads to the final form of the QCD sum rules, namely 

^\F B \ 2 e- k ' T = p{u)e-^ T dw. (11) 
I Jo 

The Borel transformation also cancels the subtraction term. The constant ratio T is called 
Borel parameter. It is an unphysical quantity in units of an energy, and the obtained values 
should be mostly independent on this parameter. This is the main criterion in analyzing 
the sum rules. 



4 Diagonal and non-diagonal correlators 

I don't want to go into details concerning the calculations of first order radiative corrections 
to the two-point correlators (for details see 0, [§). We calculate radiative corrections only 
to the leading terms in the operator product expansion because the non-leading contribu- 
tions are small. It should be stressed that we could make use of the same MATHEMATICA 
package which was constructed for the calculation of two-loop corrections of the currents it- 
self. For the case of the diagonal correlator (for an overview over the contributing diagrams 
see Fig. 2) we ended up with a renormalized spectral density 



Po orn M 



where 



with 



Po orn H 



20tH 



and r(u/fi) = t*i In ( ) + r 2 



ri = -(n 2 - 4n 7 + 6) and 
which results in the sum rule 



-F 2 (/x)e 



-A/T 



7P 



T 6 (f 5 (xc) + 



47T 



45 



(38n 2 - 137n 7 + 273 + 60C(2)), 



In 



2T 



h{xc) ~ fl{xc) ) n + f 5 (x c )r 2 



+cE 4 G T 2 f 1 (x c ) + E 6 Q exp 



2El 

J>2 



(12) 

(13) 
(14) 

+ 

(15) 



where c is a Clebsch-Gordan type factor which takes the values c = 1 for the Ag-type and 
c = —1/3 for the SQ-type doublet {Eq,Eq} ground state baryons. In order to simplify 
the notation I have introduced the abbreviations 



Ec 
T ' 



4 ' 



(Eg) 



3 ._ 



7T 



■(qq) and (E G ) 



na s (G 2 ) 



2N C ^' v 32A C (A C -1) 

The sum rule is expressed in terms of the modified Gamma functions 



(16) 



fn{x) 



X 



o ni 



c dx 



E 

rn=0 



x- 



71! 



lnx'e x dx'. 



(17) 
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I have omitted the index for the current. Later on we will see that the sum rule analysis 
is indeed independent of the chosen current within the assumed error bars. So we can use 
only one residue also for the non-diagonal case. For the non-diagonal correlator (cf. Fig. 3) 
we obtain a spectral density part 



PS ™' 



1 + ~^ r<yU l^ 



;is) 



where 



with 



pf orn (o;i 



\ u ) = — u and r y u} /^> = r i ln \2uj) 



+ r 2 



(19) 



ri = - (2n y — 8n 7 + 7 + 2{n 1 — 2)s 7 ) and 
2 

r 2 = (8n 7 -28n 7 + 37 + 8n 7 s 7 -12s 7 + 8C(2)). 



(20) 



Here the sum rule has the form 



2iVJ 



EqT 3 ( / 2 (z c ) + ( ( In ( ^ ) / 2 (rrc) - /^c) ) n + / 2 (x c )r 2 



47T 



2T 



EqEqT 1 — - fo(x c ) + - 1 — - 



c\ EqEq a s C F E, 



T 



+ 



36tt T 3 



:2r 



5 Numerical analysis 

Having the neccessary formulae at hand we can start the numerical analysis of the sum 
rules. In doing this we use the following numerical input values for the condensate contri- 
butions p. 

(qq) = -(0.23 GeV) 3 (quark condensate), (22) 
a s {G 2 ) = 0.04 GeV 4 (gluon condensate), and 
9s{(}&i±vG^ V( l) = m o{QQ) with = 0.8 GeV 2 (mixed quark-gluon condensate). 

There are in general two strategies for the numerical analysis of the QCD sum rules. The 
first strategy fixes the bound state energy A from the outset by choosing a specific value 
for the pole mass of the heavy quark and then extracts a value for the residue F. In order 
to obtain information from the sum rules which is independent of specific input values, 
we adopt a second strategy, namely to determine both A and F by finding simultaneous 
stability values for them with respect to the Borel parameter T. 

The first step in carrying out the numerical analysis of the sum rules is to find a sum 
rule "window" for the allowed values of the Borel parameter T. The parameter range of 
T is constrained by two different physical requirements. The first is that the convergence 
of the OPE expansion must be secured. We therefore demand that the subleading term 
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Baryon type state 


E c 


A 


F 


Ag (L.O.) 


1.2±0.1GeV 


0.77 ±0.05 GeV 


0.022 ± 0.001 GeV 3 


A q (N.L.O.) 


1.1 ±0.1 GeV 


0.77 ± 0.05 GeV 


0.027 ±0.002 GeV 3 


S q (L.O.) 


1.4±0.1GeV 


0.96 ±0.05 GeV 


0.031 ±0.002 GeV 3 


Eg (N.L.O.) 


1.3±0.1GeV 


0.94 ±0.05 GeV 


0.038 ±0.003 GeV 3 



Table 1: Results of the diagonal sum rule analysis for the continuum threshold parameter 
E c , the bound state energy A, and the residuum F for Ag-type and Eg-type currents, 
analyzed to leading order (L.O.) as well as next-to-leading order (N.L.O.) 

in the OPE does not contribute more than 30% of the leading order term. This gives a 
lower limit for the Borel parameter. The upper limit is determined by the requirement 
that the contributions from the excited states plus the physical continuum (even after 
Borel transformation) should not exceed the bound state contribution. This requirement 
is neccessary in order to guarantee that the sum rules are as independent as possible of the 
model-dependent assumptions concerning the profile of the theoretical spectral density, i.e. 
the model of the continuum. The lower limit of Eq is given by the requirement that the 
indicated window should be kept open. For the rest, Eq is a free floating variable which is 
only limited by the stability requirements on A and F . 

As an example I want to show the analysis of the diagonal sum rule for the Ag-type state in 
Fig. 4. Fig 4(a) shows the dependence of the bound state energy A on the Borel parameter 
T and Fig. 4(b) shows the dependence of the residue on T, both for the leading order sum 
rule. Fig. 4(c) and Fig. 4(d) show the same dependencies for the radiatively corrected sum 
rules. The same analysis is repeated for the Eg-type baryons. The results of the numerical 
analysis both without and with radiative corrections are given in Table [1] 

The numerical results for the non-diagonal sum rules are given in Table g. Assuming 
relative errors of 10% for the bound state energy and 20% for the residue, the obtained 
values are in agreement with the results of the analysis of the diagonal sum rules, where 
the values for the Eg-type baryon are the more reliable one. 

For the linear combination J = (J x ± J^jl of currents the light-side Dirac structure can 
be seen to appear in the form |(1 + i))T, i.e. one has the projector factor P + = (1 + 1)/2 
which projects on the large components of the light quark fields. We refer to this particular 
linear combination of currents as the constituent type current. This linear combination of 
currents is expected to have minimum overlap with the heavy ground state baryons in the 
constituent quark model, i.e. where the light diquark state in the heavy baryon is taken to 
be composed of on-shell light quarks. This model emerges in the large iV c -limit. 

The use of a constituent type interpolating current J = ( J\ ± J 2)12 combines the two sum 
rule formulas for the diagonal and the non-diagonal case, taking one half of each part. 
The numerical results of the analysis are given in Table |[ The constituent type sum rules 
show an improved stability on the Borel parameter T as compared to the non-diagonal 
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Baryon type state 


E c 


A 


F 


A Q (L.O.) 
Aq (N.L.O.) 

E Q (L.O.) 
E Q (N.L.O.) 


1.0 ±0.10 GeV 
1.0 ±0.10 GeV 
1.5 ±0.10 GeV 
1.2 ±0.10 GeV 


0.75 ±0.10 GeV 
0.72 ±0.10 GeV 
1.16 ±0.10 GeV 
0.94 ±0.10 GeV 


0.024 ± 0.002 GeV 3 
0.032 ±0.003 GeV 3 
0.045 ± 0.003 GeV 3 
0.039 ±0.004 GeV 3 



Table 2: Results of the non-diagonal sum rule analysis for the continuum threshold param- 
eter E c , the bound state energy A, and the residuum F for Ag-type and E^-type currents, 
analyzed to leading order (L.O.) as well as next-to-leading order (N.L.O.) 



Baryon type state 


Eq 


A 


F 


Aq (L.O.) 
A Q (N.L.O.) 

E Q (L.O.) 
E Q (N.L.O.) 


1.1 ±0.10 GeV 

1.1 ±0.10 GeV 
1.3 ±0.10 GeV 

1.2 ±0.10 GeV 


0.77 ±0.10 GeV 
0.77 ±0.10 GeV 
1.03 ±0.10 GeV 
0.94 ±0.10 GeV 


0.034 ±0.004 GeV 3 
0.032 ±0.004 GeV 3 
0.045 ±0.004 GeV 3 
0.036 ±0.004 GeV 3 



Table 3: Results of the constituent type mixed sum rule analysis for the continuum thresh- 
old parameter E c , the bound state energy A, and the residuum F for Ag-type and Eg-type 
currents, analyzed to leading order (L.O.) as well as next-to-leading order (N.L.O.) 



sum rules, but the stability is not as good as in the diagonal case. Within the assumed 
errors the results are again in agreement with both the diagonal and the non-diagonal sum 
rule analysis. 

6 Comparison with experimental values 

Taking the experimental results for charm-quark baryons, namely m(A c ) = 2284.9 ± 
0.6 MeV and m(E+) = 2453.5 ± 0.9 MeV, our central values A(A Q ) = 760 MeV and 
A(Eq) = 940 MeV predict a mean pole mass of m c = 1520 MeV for the charm quark. For 
the experimental value m(A b ) = 5642 ± 50 MeV for the mass of the baryon A b , our central 
value A(Aq) = 760 MeV for the bound state energy suggests a pole mass of m b = 4880 MeV 
for the bottom quark. 
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7 Conclusion and outlook 



• The determination of the two-loop anomalous dimension is essential for the evolution 
behaviour of the current to first order perturbation theory. 

• The presented results as well as the related computer package are applicable and 
extendable to the calculation of three-point correlators and thus the determination 
of the Isgur-Wise function. This work will be started in the near future. 

• Another extend is the calculation of correlators including massive lines within exact 
QCD. This will be done to test the convergence of the 1/toq expansion in HQET. 
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Figure Captions 



Fig. 1: Diagrams representing the first few vacuum expectation values of local op- 
erators due to the operator product expansion of the two-point correlator. 

Fig. 2: Radiative corrections to the diagonal correlator. 

(0) Lowest order two-loop contribution, 

(1) -(4) 0(a s ) three-loop contributions. 

Fig. 3: Radiative corrections to the non-diagonal correlator given by the dimension 
three condensate contribution. 

(0) Lowest order one-loop contribution, 

(1) -(4) 0(a s ) two- loop contributions. 

Fig. 4: Sum rule results on the non-perturbative parameters of the Aq as functions 
of the Borel parameter T. Shown are five curves for five different values 
of the threshold energy E c spaced by 100 MeV around the central value 
Eq = E c cst . Eq grows from bottom to top. These are in detail 

(a) lowest order sum rule results for the bound state energy A (A); 

(b) lowest order sum rule results for the absolute value of the residue F A ; 

(c) 0(a s ) sum rule results for the bound state energy A (A); 

(d) 0(a s ) sum rule results for the absolute value of the residue Fa. 
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